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Abstract - The in-plane static screening of the field originated by a charge 
placed in a graphene sheet is investigated. A self-consistent field equation in 
the real space domain is obtained by using a suitable Thomas-Fermi procedure. 
Exact and approximated (for qualitative considerations) solutions are presented. 
In the case of a charged sheet, the screened potential presents a tail dependent on 
the free carrier density whose importance is connected with the local features 
of the impurity system. Early conclusions about Thomas-Fermi screening in 
graphene are revised. 

PACS: 71.10.Ca, 73.61.Le 

Graphene is a bidimensional structure which has newly attracted the in- 
terest of physic community since the recent report on the successful isolation 
of free-standing carbon monolayers [1]. The striking properties of this mate- 
rial ensue from the peculiar electronic structure which in the ideal honeycomb 
lattice shows gapless (approximately) linear-spectrum bands with degeneracy 
points (two per cell) at the corners of the hexagonal Brillouin zone [2]. Due to 
a peculiar suppression of localization effects, electrical conductivity never falls 
below a certain value (resistivity p max « 6.5 fcO) corresponding to the quan- 
tum unit of Mott conductance [1]. Anomalous effects in transport properties 
were successfully explained in the framework of the quantum electrodynamics 
[3] which idealize the electron spectrum (near the two degeneracy points) as the 
one of a 2D-gas of massless two-fermion-Dirac quasiparticles [4]. Carrier pop- 
ulations can be changed by thermal excitation or/and by electrically induced 
band shift. Electron or holes are poured into bands in order to match the equi- 
librium Fermi level. Accordingly, electron or hole conductivities can be settled 
by applying suitable gate voltages [1]. 

This paper is devoted to the investigation of the static screening properties 
in the gapless bidimensional semiconductor. Really, screening in bidimensional 
electron gas was already investigated by using the self-consistent-field dielectric 
formulation of Ehrenreich and Cohen (EC) [5]. According to the EC's work, 
involved in 3D electron gas, the self consistency can be attained by combining 
the single-particle Liouville and the Poisson equations [6] , thus allowing for the 
Lindhard longitudinal dielectric constant [7]. The repetition of this procedure 
in a 2D system [5] may be questioned since it lacks a proper handling of Poisson 
equation for the confined electron gas. A different approach, still based on 
the EC's work, was addressed to investigate static screening in graphite layers 
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[8]. To attain the self-consistency, the static solution of the Liouville equation 
was inserted in the calculus of the induced potential. Layers was treated as 
bidimensional electron systems each having a disk-shaped Fermi (parabolic- 
spectrum) distribution with radius fci? . Both the Thomas- Fermi (TF) screening 
and the Friedel oscillations were calculated. The authors concluded that, due to 
the electron confinement, the TF screening in the layer is independent of fcp and 
that it is asymptotically ruled by a law such as V(x) ~ A cxp(—Bx)/x where A 
and B are suitable constants. The case of vanishingly small carrier density was 
handled in the short-range limit by assuming that the (RPA) screening come 
out from all the it electron (one electron per atom). Besides the latter model 
assumption, we believe that conclusions about features of the TF screening , 
that is the independence of the carrier density are to be reconsidered. 

Differently from early works, the TF screening in the graphene sheet will be 
investigated in the real space domain. To properly account for the background 
polarization the screening problem will be addressed on the line of the Oliva's 
model [9] . This model extends the TF theory to include semiconductors or insu- 
lators in which a perturbing field succeeds in "piercing" the electron forbidden 
gap, thus allowing valence-band electrons to penetrate the conduction band. 
Accordingly, the local excess of electron density is determined by the band-shift 
across the higher occupied energy level. To avoid unnecessary complications 
we disregard effects due to polarization of the medium surrounding. Of course, 
the use of the TF approach leaves out effects such as Friedel oscillations which 
become overwhelming at large distance from the source of the perturbing field. 
We will consider both the cases of a charged (by gate voltage) and uncharged 
graphene sheet. The screening problem leads to a non- linear intcgro-differcntial 
equation, but essential features of the screened potential can be deduced from 
a simplified form of such an equation. It will be shown that the decay of the 
potential depends on the carrier density in the graphene sheet. 

Let us consider a large (virtually infinite) carbon monolayer at the temper- 
ature T = K in which a gate voltage produces an electron-level energy shift 
— eU < (— e < 0). The layer is thus charged by an uniform electron distribution 
with density n* . If a point-like positive charge particle if placed in such a bidi- 
mensional system, the electron distribution adjust itself to reach the minimum 
energy state with a density n(r) dependent on the distance from the external 
charge. The induced change of charge density is thus — eSn(r) = — e (n(r) — n*). 
For simplicity's sake we will consider the bidimensional system as a continuum 
with circular isotropic properties. The dipolc background polarization effect is 
phcnomenologically subsumed in a constant e. We will return later on the afore 
presented model assumptions. Now, if we fix the the origin of the polar coor- 
dinate system on the external charge, the induced potential can be calculated 
as 
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For r > we have 
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Vmd (r) — — j p(x)Sn(xr) dx (1) 
o 

where 

p(x) = J x(x 2 -2xcos$+l)- 1/2 d$ (2) 

which shows p(0) = (dp(x = 0)/dx = 2tt), p(oo) — 2ir (dp(x = oo)/dx = 0) 
and a singularity at x=l. Since 5n(r) is absolutely integrable, the integral 
function |T]) is continuous and its derivatives can be obtained by differentiating 
under the integral sign. To avoid ambiguities in handling singular functions, 
we replace the function p{x) with a function p(x) which differ from the former 
only in a suitably small interval around x=l where it is finite and continuous as 
well as its derivatives. It is not necessary to specify the exact form of p{x) , we 
only require that — (l/27r) / °° (d 2 p(x)/dx 2 ) dx — 1 compatibly with the total 
charge-screening condition, that is, 2n Sn(r) r dr = 1. On the other hand, 
we should remove the singularity to avoid considering divergent self-interactions. 
By defining V (r) = e/r + Vmd (r) it is easy to verify that (r > 0). 
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We will simplify this equation later. Now, we search for the self-consistency by 
relating 8n(r) to V(r) , that is [9], 



5n{r) = 279 ' eM(e-eV(r)-eU)/kT] + l d£ (4) 



~ 2lg ' exp[(e-eU)/kT] + l d£ 
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where we used the density of states N(e) = 4 jer j /37T7ga 2 (70 ~ 1 eV, a = 2A6A) 
[3] and 7 9 = 4/3 7r 7ga 2 w 0.07 eV~ 2 A~ 2 The infinite integration range is 
allowed by the small energies involved in the problem dealt with. Really, we are 
also assuming that the electronic energy spectrum maintains its linear law over 
a large range of energies [1,3]. At T = we obtain 

Sn(r) =e 2 7s {[F(r) + f/] 2 -C/ 2 } (5) 

In the following, we will present the screened potential in the more convenient 
form 

V(r) = e ^ (6) 
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where the screening factor satisfies the conditions /(0) = 1/e and /(oo) = 0. 
Thus, after substitution of cq. ([5]) into eq. ([3]) and integration by parts we 
obtain 
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where the function ip(x) = — (l/2ir)xd 2 p(x)/dx 2 is everywhere small (y(0) = 
0, y(oo) = 0) and negative besides a close interval around x = 1 where it oscil- 
lates and shows its maximum positive value at x = 1. Due to these properties, in 
some favorable cases (very small U) the solution of eq. does not greatly de- 
viates from the one obtained by replacing tp(x) with the Dirac function <5(x— 1). 
This allows us to investigate qualitatively some properties of the screening fac- 
tor which will be otherwise hidden by the numerical calculus. Thus, only for 
qualitative discussion, we consider 

Equation |8|) can be analytically solved in the linear range, that is, in the range 
of distances where \V (r)| << \U\ holds so that we can write 

The solution of the linearized equation is [10] 

/ (r) = AV(4Ar)^i (V4Ar) (10) 

where K i (x) stands for the modified Bessel (Basset) function of the second kind, 
A = 4ne 3 j g U/e — (ei7/7o) 31/oe is the inverse of the "screening length" and A 
is a factor merging the short and long range solutions. We should take care in 
defining the screening length since it may happen that when linearization holds 
the potential is significantly damped by the short range screening. Really, when 
linearization holds we can use the asymptotic form of the Bessel function, that 
is [11], 

/ (r) = A (4Ar) 1/4 exp(-\/4Ar) . (11) 

In the short range distances where \V (r)| >> \U\ the form of f(r) is indepen- 
dent of U, that is, independent of free carrier density. For language convenience, 
we will refer to this case as the intrinsic screening. In the case U = , eq.. © 
becomes 



d 2 f(r) 
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where a = 2ne 4r y g /e ~ 91/e. We are able to give an asymptotic solution of this 
equation. Indeed, for f(r) << 1 we can write 
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which, when higher order term can be disregarded, has the solution 

f(r) = l/a\n(r/b) (14) 

where b is a suitable length constant whose value depend on the dielectric con- 
stant. From eq. I]14[) it is clear that screening becomes more efficient as e 
decreases. If we calculate the screening factor for two different dielectric con- 
stant ei and 62 with ti > ex, at (very large) distances r where the asymptotic 
form holds, we have /i (r)//2 (r) ~ £i/e2- 

Equation (fl"2|) is clearly scale-factor free since a is dimensionless. This spe- 
cial property cannot be ascribed only to the bidimensional character of the 
problem dealt with. But, it also come out from the linearity of the electronic 
spectrum. As a mere speculation, an analogous equation would play in the case 
of electroweak screening by the neutrino sea or, in general, by 3D gapless linear- 
spectrum fermion systems (the screening factor will appear as /(r) 3 ). Unfortu- 
nately, this property makes the screening factor sensitive of the impurity-system 
features. For more clarity, let us assume that the external charge has a radius 
u . We change the boundary conditions as f(u) — 1/e. Thus, the potential in 
the intrinsic screening case can be written also as 

VM = e t& (15) 
ux 

where x = r/u > 1. From comparison between the potential (fTS"]) and and the 
level shift induced in the charged graphene sheet we can estimate the reduced 
distance x* where the carrier density effect becomes significant, that is, x* ~ 
ef(x*)/Uu . Thus, if u — > we have x* — > oo and f(x*) — > (vanishingly 
small A values). In few words, the importance of free carrier density decreases 
as u decreases. In this connection, it is to be pointed out that several model 
assumptions fails as u — > 0, mainly: 1) the screening is not bidimensional near 
the external charge; 2) the electronic spectrum is not linear over the energies 
spanned by the large V{r); 3) the system cannot be considered as a continuum. 
However, these are general troubles and we can reasonably assume that the 
impurity as a radius not too small, namely u = 1A. The results are not 
qualitatively affected if we use the half or the twice of this value. 

For the numerical integration of eq. ([7]) we used the constraints on p(x) or, 
equivalently, the condition ip(x)dx = 0. Figure 1 shows the results for the 
cases: e = 10, U = 5 1(T 4 V (solid curve A), U = 5 1CT 6 V (solid curve B) and 
U = (solid curve B). We used small U values to facilitate numerical calculus. 
The dashed curves represent the corresponding exact solutions of eq. ©• The 
asymptotic expressions (jlll) and ([T4")) (A = 0.91 and A = 0.033 for A and B , 
respectively, and b = 7.54 A) are shown by means of dotted curves. It appears 
that deviations between solutions and ([8]) increase as U increases [13]. 

As a final note, we briefly consider the temperature effect on the screening 
charge density the case |eV(r)| << |eZ7| << kT. Now, we must take into 
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account that the induced charge density —e5n(r) includes the two contributions 
—e5ri- (r) and —eSn + (r) arising from electrons and holes, respectively. It is easy 
to verify that eq. (g]) satisfies Sn(r) = [8n-(r) — 6n+(r)] . Thus, by expanding 
the integral with respect to (eV(r) + eU)/kT and by substituting in eq. (O , 
we obtain 

Sn(r) sa 2ej g (£kT + nell) V (r) (16) 

where, by putting g(x) — \x\ exp(x)/ (exp(x) + l) 2 , £ = f^°° g(x) dx and k = 

g(x) (2exp(x) — 1) / (cxp(a;) + l)dx. Note that In the case {7=0 the 
screening parameter increases proportionally to the temperature, that is, A(T) = 
4ire 2 j g t;kT/e. The exact numerical calculus of the screening factor appear at 
this level an unnecessary complication. Thus, we do not dwell upon this point.. 



References 

a masalis@unica.it 

[I] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katshelson, 
I. V. Grigorieva, S. V. Dubonos and A. A. Firsov, Nature 438, (2005) 197; K. 
S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhnag, I. V. Grigorieva, 
S. V. Dubonos and A. A. Firsov, Science 306, (2004) 666. 

[2] P. R. Wallace, Phys. Rev. B 71 (1947) 622. 

[3] N. M. Peres, F. Guinea, A. H. Castro Neto, Phys. Rev B 73, 125411 
(2006); V. P. Gusynin and S. G. Sharapov, Phys. Rev. B 73 (2006) 245411; 

[4] G. W. Semenoff, Phys. Rev. Lett 24 (1984) 2449. 

[5] F. Stern, Phys. Rev. Lett. , 18 (1967) 546. 

[6] H Eherenreich and M. H. Cohen, Phys. Rev. 115 (1959) 786. 

[7] G. D. Mahan, Many Particle Physics (Plenum Press, New York-London) 
1981, Ch. 5 

[8] P. B. Visscher and L. M. Falicov, Phys. Rev. B 3 (1971) 2541. 
[9] J. Oliva, Phys. Rev. B 35 (1987) 3431. 

[10] I. S. Gradshteyn and I. M. Ryzhik, Tables of integrals, Series and Prod- 
ucts, Sixth Edition (Academic Press) 2000, p. 922 n.5 

[II] E. T. Wittaker and G. N. Watson, A Course of Modern Analysis (Cam- 
bridge at the University Press) 1952, p. 374 

[13] When writing this article the author was not aware of a similar paper 
published in november 2006: M. I. Katnelson, Phys. Rev. B 74, 201401 (2006). 
By applying the Fourier transforms to a linearized integral self-consistent-field 
equation it was found that f(r) decreases asymptotically as ~ 1/r 2 which, 
substantially, confirms the result of ref. [5] . 



6 



Caption of Fig. 1 

Screening factor dependence on the distance from the external charge with 
radius u = 1A in a sheet with e = 10. Solid curves represent the exact numerical 
solutions of 0: U = 5 1CT 4 V (A), U = 5 1CT 6 V (B) and U = (C). Dashed 
curves represent the exact numerical solutions of eq. |J8]): U — 5 1CP 4 V (A), 
ell = 5 1CP 6 (B) and ?7 = (C). Dotted curves represent the asymptotic 
solution of eq. © with A = 0.91 (A), A = 0.033 (B) and b = 7.54 A. (C). 
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